Measurements of CP Violation in B° — > D* iz + and B° — > D ir + Decays 

F. J 



Ronga, 7 T. R. Sarangi, 7 K. Abe, 7 I. Adachi, 7 H. Aihara, 45 D. Anipko, 1 K. Arinstein, 1 Y. Asano, 49 

T. Aushev, 11 A. M. Bakich, 40 V. Balagura, 11 E. Barberio, 19 M. Barbero, 5 K. Belous, 10 U. Bitenc, 12 I. Bizjak, 12 

S. Blyth, 22 A. Bondar, 1 A. Bozek, 25 M. Bracko, 7 ' 18 ' 12 T. E. Browdcr, 5 P. Chang, 24 A. Chen, 22 W. T. Chen, 22 

B. G. Cheon, 3 R. Chistov, 11 Y. Choi, 39 A. Chuvikov, 33 J. Dalseno, 19 M. Danilov, 11 M. Dash, 50 A. Drutskoy, 4 

S. Eidelman, 1 D. Epifanov, 1 S. Fratina, 12 N. Gabyshev, 1 A. Garmash, 33 T. Gershon, 7 G. Gokhroo, 41 B. Golob, 17 ' 12 

A. Gorisek, 12 J. Haba, 7 K. Hara, 7 T. Hara, 30 K. Hayasaka, 20 H. Hayashii, 21 M. Hazumi, 7 T. Higuchi, 7 L. Hinz, 16 

T. Hokuue, 20 Y. Hoshi, 43 S. Hou, 22 W.-S. Hou, 24 T. Iijima, 20 K. Inami, 20 A. Ishikawa, 45 H. Ishino, 46 R. Itoh, 7 

Y. Iwasaki, 7 J. H. Kang, 51 H. Kawai, 2 T. Kawasaki, 27 H. R. Khan, 46 H. J. Kim, 15 K. Kinoshita, 4 S. Korpar, 18 ' 12 

P. Krokovny, 1 R. Kulasiri, 4 R. Kumar, 31 C. C. Kuo, 22 Y.-J. Kwon, 51 G. Leder, 9 J. Lee, 37 T. Lesiak, 25 J. Li, 36 

A. Limosani, 7 D. Liventsev, 11 G. Majumder, 41 F. Mandl, 9 T. Matsumoto, 47 W. Mitaroff, 9 K. Miyabayashi, 21 

H. Miyake, 30 H. Miyata, 27 D. Mohapatra, 50 T. Nagamine, 44 I. Nakamura, 7 E. Nakano, 29 Z. Natkaniec, 25 

S. Nishida, 7 O. Nitoh, 48 S. Noguchi, 21 T. Nozaki, 7 S. Ogawa, 42 T. Ohshima, 20 S. Okuno, 13 S. L. Olsen, 5 

Y. Onuki, 27 P. Pakhlov, 11 C. W. Park, 39 H. Park, 15 L. S. Peak, 40 R. Pestotnik, 12 L. E. Piilonen, 50 A. Poluektov, 1 

M. Rozanska, 25 Y. Sakai, 7 N. Sato, 20 N. Satoyama, 38 K. Sayeed, 4 T. Schietinger, 16 O. Schneider, 16 A. J. Schwartz, 4 

R. Seidl, 6 - 34 K. Senyo, 20 M. E. Sevior, 19 M. Shapkin, 10 H. Shibuya, 42 B. Shwartz, 1 J. B. Singh, 31 A. Sokolov, 10 

A ^omov, 4 N. Soni, 31 R. Stamen, 7 S. Stanic, 28 M. Staric, 12 H. Stoeck, 40 K. Sumisawa, 30 S. Suzuki, 35 S. Y. Suzuki, 7 

F. Takasaki, 7 M. Tanaka, 7 Y. Tcramoto, 29 X. C. Tian, 32 K. Trabelsi, 5 T. Tsukamoto, 7 S. Uehara, 7 

K. Ueno, 24 S. Uno, 7 P. Urquijo, 19 Y. Ushiroda, 7 Y. Usov, 1 G. Varner, 5 K. E. Varvell, 40 S. Villa, 16 

C. H. Wang, 23 Y. Watanabe, 46 E. Won, 14 Q. L. Xie, 8 B. D. Yabsley, 40 A. Yamaguchi, 44 Y. Yamashita, 26 

M. Yamauchi, 7 J. Ying, 32 C. C. Zhang, 8 L. M. Zhang, 36 Z. P. Zhang, 36 V. Zhilich, 1 and D. Ziircher 16 

(The Belle Collaboration) 

1 Budker Institute of Nuclear Physics, Novosibirsk 

2 Chiba University, Chiba 

" 'Chonnam National University, Kwangju 

' University of Cincinnati, Cincinnati, Ohio 45221 

5 University of Hawaii, Honolulu, Hawaii 96822 

6 University of Illinois at Urbana- Champaign, Urbana, Illinois 61801 

7 High Energy Accelerator Research Organization (KEK), Tsukuba 

8 Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 

9 Institute of High Energy Physics, Vienna 

10 Institute of High Energy Physics, Protvino 

irptirnl nnrf J^'rnfnrTnfritnl PhysiQ§ MOSCOW 



o 
o 

(N 

Sh 

<: 

A. Somov, 4 N. Soni, 3 
F. Takasaki, 7 

> 

o 
o 

O 

*s 

<D 

43 



Y**~ riLSLiLULe oj ntgu rjiieigy rtiysics, riot 

IJ^l ' 11 Institute for Theoretical and Experimental Phy. 

rS \ 12 J. Stefan Institute, Ljubljana 

Jj . ls Kanagawa University, Yokohama 



Kanagawa University, Yokohama 
"Korea University, Seoul 

Taegu 



4 Korea University, Seoul 
J Kyungpook National University, 
'Swiss Federal Institute of Technology of Lausanne, EPFL, Lausanne 
17 University of Ljubljana, Ljubljana 

18 University of Maribor, Maribor 
19 University of Melbourne, Victoria 

20 Nagoya University, Nagoya 

21 Nara Women's University, Nara 

22 National Central University, Chung-li 

23 National United University, Miao Li 

^Department of Physics, National Taiwan University, Taipei 

25 H. Niewodniczanski Institute of Nuclear Physics, Krakow 

Nippon Dental University, Niigata 

27 Niigata University, Niigata 

28 Nova Gorica Polytechnic, Nova Gorica 

29 Osaka City University, Osaka 

Osaka University, Osaka 
31 Panjab University, Chandigarh 
Peking University, Beijing 
33 Princeton University, Princeton, New Jersey 08544 



H RIKEN BNL Research Center, Upton, New York 11973 
Saga University, Saga 

36 i 



' University of Science and Technology of China, Hefei 
87 Seoul National University, Seoul 

J8 Shinshu University, Nagano 

39 Sungkyunkwan University, Suwon 

**° University of Sydney, Sydney NSW 

^ Tata Institute of Fundamental Research, Bombay 

^ 2 Toho University, Funabashi 

^Tohoku Gakuin University, Tagajo 

™ Tohoku University, Sendai 

■* Department of Physics, University of Tokyo, Tokyo 

Jtb Tokyo Institute of Technology, Tokyo 

^ Tokyo Metropolitan University, Tokyo 

Jt8 Tokyo University of Agriculture and Technology, Tokyo 

^University of Tsukuba, Tsukuba 

50 Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24-061 

jl Yonsei University, Seoul 

We report measurements of time dependent decay rates for B° — > D'*- >_ 7r + decays and extraction 
of CP violation parameters that depend on 03. Using fully reconstructed D^*'ir events and partially 
reconstructed D*n events from a data sample that contains 386 million BB pairs that was collected 
near the T(45) resonance, with the Belle detector at the KEKB asymmetric energy e + e~ collider, 
we obtain the CP violation parameters S + (D ( * ) tt) and S~(D m tt). We obtain S + (D*tt) = 0.049 ± 
0.020(stat) ± O.Oll(sys), S~(D*tv) = 0.031 ± 0.019(stat) ± O.Oll(sys), and S + (D%) = 0.031 ± 
0.030(stat)±0.012(sys), S~(Dn) = 0.068 ± 0.029(stat) ±0.012(sys). These results are an indication 
of CP violation in B° — » D*~tt + and B° — > D~tt + decays at the 2.5 a and 2.2 a levels, respectively. If 
we use the values of R D (*)^ that are derived using assumptions of factorization and SU(3) symmetry, 
the branching fraction measurements for the D a * -k modes, and lattice QCD calculations, we can 
restrict the allowed region of | sin(20i + 03)| to be above 0.44 and 0.52 at 68% confidence level from 
the D*n and Dir modes, respectively. 



PACS numbers: 13.65.+i, 13.25.Gv, 14.40.Gx 



I. INTRODUCTION 



Within the Standard Model (SM), CP violation arises 
due to a single phase in the Cabibbo-Kobayashi-Maskawa 
(CKM) quark mixing matrix [1]. Precise measurements 
of CKM matrix parameters therefore constrain the SM, 
and may reveal new sources of CP violation. Measure- 
ments of the time-dependent decay rates of B°(B°) — > 
£)(*)T 7r ± provide a theoretically clean method for ex- 
tracting sin(20i + 03 ) [2]. As shown in Fig. 1, these 
decays can be mediated by both Cabibbo-favoured decay 
(CFD) and doubly-Cabibbo-suppressed decay (DCSD) 
amplitudes, V* b V u d and V* b V c d, which have a relative 
weak phase 03. 





(a) CFD 



(b) DCSD 



FIG. 1: Diagrams for (a) B" -> D { *>--k + and (b) B 
£) (*)- 7r +. Those for b° -+ D (,)+ tt- and B° -> D ( * )+ 7r^ 
be obtained by charge conjugation. 



The time-dependent decay rates are given by [3] 
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Here At is the difference between the time of the decay 
and the time that the flavour of the B meson is tagged, 
tbo is the average neutral B meson lifetime, Am is the 
B°-B° mixing parameter, and C = (l - R 2 ) / (l + R 2 ), 
where R is the ratio of the magnitudes of the DCSD and 
CFD (we assume the magnitudes of both the CFD and 
DCSD amplitudes are the same for B° and B° decays). 
The CP violation parameters are given by 
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where L is the orbital angular momentum of the final 
state (1 for D*ir and for Dir), and 5 is the strong phase 
difference of the CFD and DCSD. The values of R and 5 
are not necessarily the same for D*ir and Dir final states, 
and are denoted with subscripts, D*ir and Dir, in what 
follows. 

Although not measured yet, the value of R is predicted 
to be about 0.02 [4]. Therefore, we neglect terms of 
O (R 2 ) (and hence take C = 1). There are theoreti- 
cal arguments that the still unmeasured values of S for 
both D*tt and Dir are small [3, 5]. 

Due to the size of R, CP violation is expected to be 
a small effect in these decays. Therefore, a large event 
sample is needed in order to obtain sufficient sensitiv- 
ity. With this in mind, we employ a partial reconstruc- 
tion technique [6] for the D*ir analysis in addition to 
the conventional full reconstruction method. In this ap- 
proach, the signal is distinguished from background on 
the basis of kinematics of the "fast" pion from the de- 
cay B — ► D*nf, and the "slow" pion from the decay 
D* — ► Dir s , alone; no attempt is made to reconstruct the 
D meson from its decay products. Background from con- 
tinuum e + e~ — ► qq (q = u, d, s, c) events is reduced dra- 
matically by requiring the presence of a high-momentum 
lepton in the event, which also serves to tag the flavour 
of the associated B in the event. 

Results from our previous analyses using a 140 fb _ 
data sample containing 152 million BB pairs have been 
published for the full reconstruction method [7] and 
the partial reconstruction method [8]. Results from 
similar analyses by BaBar collaboration were also re- 
ported [9, 10]. This study is a continuation of similar 
analyses with a substantially increased data sample con- 
taining 386 million BB events, and several improvements 
in the analyses. 



II. BELLE DETECTOR 

The data was collected with the Belle detector [11] 
at the KEKB asymmetric energy electron-positron col- 
lider [12] operating at the T(4S) resonance. The Belle de- 
tector is a large-solid-angle magnetic spectrometer that 
consists of a silicon vertex detector (SVD), a 50-layer 
central drift chamber (CDC), an array of aerogel thresh- 
old Chcrcnkov counters (ACC), a barrel- like arrangement 
of timc-of- flight scintillation counters (TOF), and an 
electromagnetic calorimeter (ECL) comprised of CsI(Tl) 
crystals located inside a superconducting solcnoidal coil 
that provides a 1.5 T magnetic field. An iron flux-return 
located outside of the coil is instrumented to detect K\ 
mesons and to identify muons (KLM) . A sample contain- 
ing 152 million BB pairs was collected with a 2.0 cm 
radius bcampipe and a 3-layer silicon vertex detector 
(SVD1), while a sample of 234 million BB pairs was col- 
lected with a 1.5 cm radius beampipe, a 4-layer silicon 
detector, and a small-cell inner drift chamber (SVD2) 
[13]. 



III. FULL RECONSTRUCTION ANALYSIS 



A. Signal Event Selection 



The selection of hadronic events is described else- 



D* + ir event selec- 



where [14]. For the B° 
tion, we use the decay chain D* 
D°tt + with subsequent decays of D + — ► K~n + ir + and 
D° -► K--K+, R-ir+n , K^tt+tt+t:-, K s ir + ir-(K s -> 
7r + 7r~). (Charge conjugate modes are implied through- 
out this Paper.) All charged tracks except for the slow 
pions from D* — ► Dir decays are required to have a min- 
imum of one hit (two hits) in the r-<j> (z) coordinate of 
the vertex detector, where the r-<p plane is transverse to 
the positron beam line that defines the z axis. These 
requirements allow a precise determination of the pro- 
duction point. To separate kaons from pions, we form 



likelihood for each track, C 



K{-k)- 



The kaon likelihood 



ratio, P(K/ir) = Ck/{C-k + C„), has values between 
(likely to be a pion) and 1 (likely to be a kaon). We 
require charged kaons to satisfy P(K/ir) > 0.3, corre- 
sponding to about 95% efficiency for detecting kaons and 
about 2% probability for misidentifying pions as kaons. 
There is no such requirement for charged pions coming 
from D decays. 

Neutral pions are formed from photon pairs with in- 
variant masses between 0.118 GeV/c 2 and 0.150 GeV/c 2 . 
The photon momenta are then recalculated with a ir° 
mass constraint. 



We require the invariant mass of K % 



r + 7r candi- 
dates to be between 0.485 GeV/c 2 and 0.510 GeV/c 2 cor- 
responding to ±5cr, where a is the Monte Carlo (MC) 
determined invariant mass resolution. The radial impact 
parameter, which is the distance of closest approach of 
the candidate charged tracks to the event-dependent in- 
teraction point (IP) in the r-<p plane, is required to be 
larger than 0.02 cm for high momentum (> 1.5GeV/c) 
Kg candidates and 0.03 cm for those with momentum 
less than 1.5GeV/c. Here the IP is determined for each 
set of 10,000 neighboring hadronic events, and the Kg 
momentum is given in the T(4S) rest frame (cms). The 
7r + 7r~ vertex is required to be displaced from the IP by 
a minimum transverse distance of 0.22 cm for high mo- 
mentum candidates and 0.08 cm for the remaining candi- 
dates. The mismatch in the z direction at the Kg vertex 
point for the ir + 7r~ tracks must be less than 2.4 cm for 
high momentum candidates and 1.8 cm for the remain- 
ing candidates. The direction of the pion pair momen- 
tum must also agree with the direction of the vertex point 
from the IP to within 0.03 rad for high momentum candi- 
dates and to within 0.1 rad for the remaining candidates. 
A mass and vertex constraint is imposed when fitting the 
Kg candidates. 

For D° meson candidates, the invariant 
mass of the D° candidate is required to be 
within ±20, ±30, ±20, and ±20MeV/c 2 for 
K~tt + , K~ir + ir°, K~ir + ir + ir~ , K s ir + ir~ modes 
respectively, while the invariant mass of the D + candi- 



dates should be within ±18MeV/c 2 of the nominal D + 
mass corresponding to ±4cr. For the D° — » _ftr~7r + 7r° 
mode, we further require the ir° cms momentum to 
be greater than 200MeV/c. We use a mass- and 
vertex-constrained fit for D candidates. 

The D* + is reconstructed by combining either D° can- 
didates with a slow tt + meson, or D + candidates with a 
slow tt° meson. The D* candidates are required to have a 
mass difference AM — Md-k ~ Md within ±3 MeV/c 2 to 
±5 MeV/c 2 of the nominal value depending on the decay 
mode. 

We reconstruct a B candidate by combining the D^*' + 
candidate with a 7r~~ candidate satisfying P(K/ir) < 0.8, 
corresponding to more than 90% efficiency for detect- 
ing pions and less than 10% probability for misidenti- 
fying kaons as pions. We identify B decays based on 
requirements on the energy difference AE = £\ Ei — 
E hcam and the bea m-energy constrained mass Mbc = 
V^bcam ~ (12iPi) 2 > where E hcam is the beam energy, 
and pi and Ei are the momenta and energies of the 
daughters of the reconstructed B meson candidate, all 
in the cms. We define a signal region in the AE- 
M bc plane of 5.27GeV/c 2 < M bc < 5.29GeV/c 2 and 
| AE\ < 0.045 GeV, corresponding to about ±3<r on each 
quantity. If more than one B candidate is found in the 
same event, we select the one with best D vertex quality. 
For the determination of background parameters, we use 
events in a sideband region defined by M^ c > 5.2 GeV/c 2 
and -0.14 GeV < AE < 0.20 GeV, excluding the signal 
region. 

The AE and M bc distributions for D*ir and Dir candi- 
dates are shown in Fig. 2. A study using MC events indi- 
cates the presence of "peaking background" that peaks in 
the signal AE-M^ C region and amounts to 1.7% (0.7%) 
of the D*ir (Dir) candidates, respectively. We treat this 
contribution as a part of the signal, and assign a system- 
atic error to account for this. 
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FIG. 2: AE and M bc distributions for 5° -» D*~tt + can- 
didates (a and b), and B° — > D~ir + candidates (c and d). 
Curves are the fit results. Hatched regions indicate the back- 
ground components of the fits. 



B. Flavour Tagging 

Charged leptons, pions, and kaons that are not asso- 
ciated with the reconstructed D^*>tt decays are used to 
identify the flavour of the accompanying B meson. The 
algorithm [15] produces two parameters for each event, q 
and r, where q = +l(q = —1) for a B° (B°) meson and r 
is a quality factor ranging from for no flavour discrim- 
ination to 1 for unambiguous flavour assignment. The 
algorithm uses the kinematical variables of the event and 
compares to those from a large number of MC events, 
and is used only to sort the data into six intervals of r 
according to estimated flavour purity. More than 99.5% 
of the events are assigned non-zero values of r. 



C. Vertex Measurement 

The decay vertex of the B — > D^*'n candidate is fitted 
using the track information of the D and it (except the 
slow it from D* decay). For the decay vertex of the tag- 
ging B meson, the remaining well-reconstructed tracks in 
the event are used. Tracks that are consistent with K s 
decay are rejected. We impose the additional require- 
ment that both signal-side and tagging-side vertices be 
consistent with the run-dependent IP profile. A study 
using a MC event sample shows that the nominal \ 2 of 
the vertex reconstruction is strongly correlated with the 
distance from the IP to the reconstructed vertex. To 
avoid possible bias in the event selection due to this cor- 
relation, we introduce a quantity that only depends on 
the z coordinate quantities; 
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where n is the number of tracks, Zi are the z vertices of 
the i-th track before and after the vertex fit, and e|? re 
is the measurement error of the z before the vertex fit. 
This quantity is calculated for the signal- and tagging- 
side separately, only for cases with multiple tracks. We 
require £ < 100 to eliminate badly reconstructed ver- 
tices, which amount to 4% (2%) of the signal-side and 
1% (1%) of the tagging-side in the D*-e and Dn modes, 
respectively. 

The proper-time difference between the vertices z rcc 
and 2tag (measured along the beam line) of the fully re- 
constructed B candidate and the tagging B meson, re- 
spectively, is calculated as 



At = (z rcc - z tag )//37c, 



(4) 



where [fa/ = 0.425 is the Lorentz boost factor of the centre 
of mass frame at KEKB. After application of the event se- 
lection criteria and the requirement that both B mesons 
have well defined vertices and |Ai| < 70 ps (~ 45r B o), 
31491 and 31725 candidates remain in the D*tt and Dn 



modes, respectively. The signal fractions of the candi- 
dates, which vary for different r bins, are 89% for D*ir 
and 83% for Dn. 



D. At Fit 

Unbinncd maximum likelihood fits to the four time de- 
pendent decay rates are performed to extract S + and S~ . 
We minimize —2 ^ In Li where the likelihood for the i-th 
event is given by 

A = (1 " /ol) [fsigPsig + (1 - /sig)Pbkg] + /ol-Pol- (5) 

The signal fraction / s ; g is determined from the (AE-M^c) 
value of each event. The signal distribution is the product 
of the sum of two Gaussians in AE and a Gaussian in 
Mb c ; that for the background is the product of a first 
order polynomial in AE and an ARGUS function [16] in 
M bc . 

The signal At distribution is given by 
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for the g = — 1 sample, and 
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for the <? = +1 sample. Here w + and w_ are respectively 
the probabilities to incorrectly assign the flavour of tag- 
ging B° and B° mesons (wrong tag fractions), and the 
decay rates are given by Eq. 1. 

The corresponding background distribution is param- 
eterized as a sum of a J-function component and an ex- 
ponential component with lifetime Tbk g : 



5(At - f&J + 



I 1 /bke) -\At-uZ 



2Tbka 



bW |At-^ kg |/r bkg 



(8) 



where /^ k is the fraction of events contained in the 5- 

function, ^b k g and /x^kg are the mean values of At in 
the 5-function and exponential components, respectively 
These parameters are determined separately from a fit to 
the At distribution in the AP-Mb c sideband data for the 
D*ir and Dir data samples, SVD1 and SVD2 data, and 
cases where the vertices are reconstructed using single 
track and multiple tracks. Typically, values of the back- 
ground At parameters are: /^ k w 0.4, Tbkg ~ 1.1 ps, 

A^bkg ~ -0- 1 P s > and Mbkg ~ °- 

A small number of events have poorly reconstructed 
vertices resulting in very broad At distributions. We 
account for this outlier contribution by adding a Gaus- 
sian component P \ with a width and fraction determined 
from the B lifetime analysis [23]. 



E. At Resolution 

The probability density functions (PDF) for the signal 
and background must be convolved with corresponding 
At resolution functions, that are determined on an event- 
by-event basis using the estimated uncertainties on the 
z vertex positions [17, 23], in order to be compared with 
the data. The signal resolution function is a convolution 
of three contributions: resolution functions for vertex re- 
construction, smearing due to non-primary tracks (Kg 
and charm daughters), and the kinematic approximation 
that the B mesons are at rest in the cms. The resolution 
function is described in detail elsewhere [23]. 

The background resolution function is parametrized as 
a sum of two Gaussians where different values are used 
for the parameters depending on whether or not both 
vertices are reconstructed with multiple tracks. These 
parameters are determined from the AP-Mb c sideband 
data. 



F. Tagging Side CP Violation Effect 

While the tagging side should have no asymmetry if 
the flavour is tagged by primary leptons since semilep- 
tonic decays are flavour-specific processes, it is possible 
to introduce a small asymmetry when daughter particles 
of hadronic decays such as D^*'n are used for the flavour 
tagging, due to the same CP violating effect that is the 
subject of this paper [18]. This effect is taken into ac- 
count by replacing the S parameters in Eq. 1 by 

B o ^ D (*)+ n - 

B a -> D^-n+ 
B a -> D«+7T- 
W -> D«-7T + : S- -> (S- - 5 t - ag ), (9) 

respectively. Here S^„ and «S t ~ represent the CP viola- 
tion effect on the flavour tagging side due to the pres- 
ence of B° — > D X and B° — > DX amplitudes, re- 
spectively. Note that unlike the S parameters, which 
are defined rigorously in terms of B° — ► D^'^tt and 
B° — ► D^^n^ amplitudes, S^ g are effective quantities 
that include effects of the fraction of B — > DX(D X) 
components in the tagging B decays and the subsequent 
behaviour of D(D) mesons. Therefore, these quantities 
must be determined experimentally. 

The values of S ta „ are determined in each r bin by 
fitting the At distributions of a B — > D*lv control sam- 
ple [15] using the signal PDFs of Eqs. 6 and 7 and set- 
ting S^ to zero. Since the D*lv final states have spe- 
cific flavour, any observable asymmetry must originate 
from the tagging side. The results for each r bin are 
listed in Table I. The errors listed here are statistical 
only. The result for the combined r bins arc S 1 ^ = 
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-> (6 — «3 tag ), 


s~ - 


-> (b + o tag ), 


s+- 


-► (S + O t ag) i 


s- - 
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-0.002 ± 0.009 ± 0.006 and SZ 



tag 

0.017 ±0.009 ±0.006, 



where the second errors are systematic. 



TABLE I: Parameters describing tagging side CP violation 
effect that are determined from the D*lv data sample. 
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G. Fit Result 



FIG. 3: (a) At distribution for the signal candidates 
when all four signal categories are combined. (b) 

(OF - SF)/(OF + SF) asymmetries. Curves are the fit re- 
sults. The hatched region in (a) indicates the background. 



The procedures for At determination and flavour tag- 
ging are tested by extracting rgo and Am. When all 
four signal categories in Eq. 1 are combined, the signal 
At distribution reduces to an exponential lifetime distri- 
bution as shown in Fig. 3(a). We do a simultaneous fit to 
the SVD1 and SVD2 samples by combining the D*ir and 
Dtt candidate events and obtain t^o = 1.532 ± 0.013 ps, 
where the error is statistical only, in good agreement 
with the world average of 1.536 ± 0.014 ps [19]. Com- 
bining the two CFD-dominant modes (denoted as OF 
because the signal-side and tagging-side have opposite B 
flavour) and the two mixing-dominant modes (denoted as 
SF because the signal-side and tagging-side have same B 
flavour) and ignoring the CP violating terms, an asym- 
metry (OF - SF)/(OF + SF) behaves as cos(AmAt) as 
shown in Fig. 3(b). A similar fit to the one used for the 
lifetime determination gives Am — 0.494 ± 0.007 ps -1 , 
where the error is statistical only, also in good agreement 
with the world average 0.502 ± 0.007 ps" 1 [19]. This fit 
is also used to determine the signal At resolution param- 
eters and wrong tag fractions iu+ and w- in each r bin 
for both SVD1 and SVD2 samples. Table If lists the fit 
result for the wrong tag fractions in the SVD2 sample. 
Those from the SVD1 sample have very similar values. 



TABLE II: Fit results for the wrong tag fractions of the SVD2 
sample. They are determined in each r bin as the averages 
and differences of the two flavours. 



(w+ + w-)/2 (w + 



0.000 - 0.250 
0.250 - 0.500 
0.500 - 0.625 
0.625 - 0.750 
0.750 - 0.875 
0.875 - 1.000 



0.461+0.006 
0.327+0.009 
0.216+0.010 
0.163+0.009 
0.124+0.009 
0.032+0.005 



+0.001 ± 0.009 
-0.026 + 0.013 
+0.032 + 0.015 
-0.001 + 0.014 
-0.020 + 0.014 
+0.012 + 0.009 



We then perform fits to determine the values of S 
by fixing r^o and Amj to the world average values and 
using previously determined w± and S tag in each r bin. 
The results are S+(D*n) = 0.050 + 0.029, S~(D*n) = 
0.028 ± 0.028, S+(Dn) = 0.031 ± 0.030, and S~ {D-k) = 
0.068 ± 0.029. The errors are statistical only. 

The At distributions for the subsamples having the 



best quality flavour tagging (0.875 < r < 1.0) are shown 
in Fig. 4 for the D*ir mode and in Fig. 5 for the Dtt 
mode. For both modes, this region constitutes roughly 
14% of the total, but has significant statistical power for 
the S determination. 



H. Systematic Error 

The systematic errors come from the uncertainties of 
parameters that are constrained in the fit (At resolution, 
background At shape, background fractions, wrong tag 
fractions, vertexing, physics parameters), uncertainties 
due to the tag side asymmetry and biases induced by 
the fitting method. To estimate contributions from the 
fit parameter uncertainties, we repeat the fit by varying 
each parameter by a given amount, and assign the shift in 
the S parameters from the nominal fit as a systematic 
error. The signal At resolution parameters are varied by 
±1 a of their errors. The background At shape param- 
eters are varied by ±1 a of their errors. For the back- 
ground fraction, in addition to varying the parameters of 
the AE-M^c signal region fit by ±1 a, we vary the AE 
signal region cut by ±5 MeV and M^ c signal region cut 
by ±3 MeV/c 2 and the quadratic sum of these is assigned 
as an error. Contributions from the peaking background 
are found to be negligible based on a MC study. We vary 
the wrong tag fraction parameters by ±1 a in each r bin 
and add in quadrature. We vary the cut of £ by +100 
and -50, which gives 0.004 for D*n and 0.002 for Dtt 
as errors due to the vertexing. We repeat the fit after 
removing the |At| < 70 ps cut and find a negligible shift. 
We vary the errors for tb° and Am by ±1 a. 

For the tag side asymmetry, the quadratic sum of sta- 
tistical and systematic errors of S^ parameters are var- 
ied by ±1 <j in each r bin and the deviations are added 
in quadrature, giving 0.005 as an error. Since this error 
includes contributions from unknown effects in the vertex 
measurement such as the misalignment of the SVD and 
drift of the IP profile, we do not explicitly assign addi- 
tional error to the vertexing other than that from the £ 
cut. 

Fit bias is tested using a D*ir signal MC sample where 
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FIG. 4: At distributions for the D*tt events in the 0.875 < 
r < 1.0 flavour tagging quality bin. (a) B° — > D* + tt~, (b) 
5° -> D*~-k+, (c) ~B° -» D* + ^r-, (d) B° -> D*"tt + . Curves 
show the results of fits to the entire event sample, and hatched 
regions indicate the backgrounds. 
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FIG. 5: At distributions for the Dn events in the 0.875 < 
r < 1.0 flavour tagging quality bin. (a) B° — » D + iv~ , (b) 
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show the results of fits to the entire event sample, and hatched 
regions indicate the backgrounds. 



one B meson decays to D* ir, and the other B decays 
generically. We fit the S^ parameter to this sample, with 
and without the S tag correction given in Eq. 9, taking 
S^ values from aB^ D*lv signal MC sample. Neither 
MC sample includes the tag side CP violation effects, so 
the two fits should return the same S values in princi- 
ple. We take the difference between the fits, 0.010, as a 
measure of biases in the fitting procedure. 

We obtain a total systematic error of 0.013 for D*ir and 
0.012 for Dn. Table III summarizes the contributions to 
the systematic errors. 



TABLE III: Summary of the systematic errors in the S mea- 
surements using the full reconstruction method. 



Sources 


D*TY 


Dty 


Signal At resolution 


0.005 


0.004 


Background At shape 


negligible 


negligible 


Background fraction 


0.002 


0.001 


Wrong tag fraction 


0.002 


0.002 


Vertexing 


0.004 


0.002 


Physics parameters (Am 


T B o) 0.001 


0.001 


Tag side asymmetry 


0.005 


0.005 


Fit bias 


0.010 


0.010 


Total 


0.013 


0.012 



I. Result 



We obtain 



S + (D*n) = 0.050 ±0.029 ±0.013, 

S-{D*ir) = 0.028 ±0.028 ±0.013, 

S + (Dtt) = 0.031 ±0.030 ±0.012, 

S~{Dtt) = 0.068 ±0.029 ±0.012 



(10) 



from the full reconstruction method, where the first error 
is statistical and the second error is systematic. 

IV. PARTIAL RECONSTRUCTION ANALYSIS 

A. Signal Event Selection 

Candidate events are selected by requiring the presence 
of fast pion and slow pion candidates. In order to obtain 
accurate vertex position determinations, fast pion candi- 
dates are required to have a radial (longitudinal) impact 
parameter dr < 0.1cm (\dz\ < 2.0cm), to have associ- 
ated hits in the SVD, and to have a polar angle in the lab- 
oratory frame in the range 30° < 6\ a b < 135°. The vertex 
positions are obtained by fits of the candidate tracks with 
the IP. Fast pion candidates are required to be incon- 
sistent with the lepton hypothesis (see below), and the 
kaon hypothesis, based on information from the CDC, 
TOF and ACC. A requirement on the fast pion cms mo- 
mentum of 1.83 GeV/c < p nf < 2.43 GeV/c is made; this 
range includes both signal and sideband regions (defined 
below). Slow pion candidates are required to have cms 
momentum in the range 0.05 GeV/c < p„ 3 < 0.30 GeV/c. 
No requirement is made on particle identification for slow 
pions; since they are not used for vertexing, only a loose 
requirement that they originate from the IP is made. The 
fast and slow pion candidates must have opposite charges. 

B. Flavour Tagging 

In order to tag the flavour of the associated B me- 
son and to reduce background from continuum e + e~ — > 



qq (q = u,d,s,c) processes, we require the presence of 
a high- momentum lepton in the event. Tagging lepton 
candidates are required to be positively identified either 
as electrons, on the basis of information from the CDC, 
ECL and ACC, or as muons, on the basis of information 
from the CDC and the KLM. They are required to have 
momentum in the range 1.2 GcV/c < pi ta < 2.3 GeV/c, 
and to have an angle with the fast pion candidate that 
satisfies —0.75 < cos S nf i in the cms. The lower bound on 
the momentum and the requirement on the angle also re- 
duce, to a negligible level, the contribution of leptons pro- 
duced from semi-lcptonic decays of the unreconstructed 
D mesons in the B° — > D* =F 7r ± decay chain. No other 
tagging lepton candidate with momentum greater than 
1.0 GeV/c is allowed in the event to reduce the mistag- 
ging probability, and also to reduce the contribution from 
leptonic charmonium decays. Identical vertexing require- 
ments to those for fast pion candidates are made in order 
to obtain an accurate Ztag position. To further suppress 
the small remaining continuum background, we impose 
a loose requirement on the ratio of the second to zeroth 
Fox- Wolfram [20] moments, i?2 < 0.6. 



C. Kinematic Fit 

Signal events are distinguished from background using 
three kinematic variables, which are approximately inde- 
pendent for signal. These are denoted by p 7r/ , cosd Wf7T3 
and cos^hci- For signal, the fast pion cms momentum, 
Pit-f , has a uniform distribution, smeared by the experi- 
mental resolution, as the fast pion is monoenergetic in the 
B rest frame. The cosine of the angle between the fast 
pion direction and the opposite of the slow pion direction 
in the cms, cos<5„- /7rs , peaks sharply at +1 for signal, as 
the slow pion follows the D* direction due to the small 
energy released in the D* decay. The cosine of the angle 
between the slow pion direction and the opposite of the 
B direction in the D* rest frame, cos^hei, has a distribu- 
tion proportional to cos8^ cl for signal events, as the B 
decay is a pseudoscalar to pseudoscalar vector transition. 
Since the D* is not fully reconstructed, cos #hci is calcu- 
lated using kinematic constraints, and the background 
can populate the unphysical region |cos6*h i| > 1- 

We select candidates that satisfy 1.93 GeV/c < p 7r/ < 
2.43 GeV/c, 0.850 < cosS 7Tf7Ts < 1.000 and -1.70 < 
cos^hd < 1-80. In the cases where more than one can- 
didate satisfies these criteria, we select the one with the 
largest value of cos S 7Tf7Ts . We further define signal regions 
in p^ and cosS- Rf T Ts as 2.13 GcV/c < p 7Tf < 2 A3 GeV/c, 
0.925 < cos6 nfne < 1.000, and two regions in cos6>hci: 
-1.00 < cos6» hcl < -0.30 and +0.40 < cos6» hcl < +1.10. 

Background events are separated into three categories: 
D* T p ± , which is kinematically similar to the signal; cor- 
related background, in which the slow pion originates 
from the decay of a D* that in turn originates from the 
decay of the same B candidate as the fast pion candi- 
date (e.g., D**tt); and uncorrelated background, which 



includes everything else (e.g., continuum processes, Dtt). 
The kinematic distributions of the background categories 
and the signal are determined from a large MC sample, 
corresponding to two times the integrated luminosity of 
our data sample, in which the branching fractions of the 
signal and major background sources are weighted ac- 
cording to the most recent knowledge [19, 21]. We also 
use this MC sample for various tests of the analysis al- 
gorithms. 

Event-by-event signal fractions are determined from 
binned maximum likelihood fits to the three-dimensional 
kinematic distributions (6 bins of p n . x 5 bins of cos 5„.„ a 
x 10 bins of cos^hci). The results of these fits, projected 
onto each of the three variables, are shown in Fig. 6, and 
summarized in Table IV. 
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FIG. 6: Results of the kinematic fits to D*n candidates, pro- 
jected onto (top) thepTr, axis, (middle) the cosJ^,^ s axis and 
the (bottom) cos 8hci axis, in both (left) selection region and 
(right) signal region. The arrows show the edges of the signal 
regions. 



TABLE IV: Summary of the results of the three-dimensional 
fits to kinematic variables. The numbers of events and frac- 
tions given for each category are those extrapolated to inside 
the signal regions in all three variables. The errors on the 
fractions are returned by the fit and propagated to the corre- 
sponding numbers of events. 

Candidates Fraction 



Data 

D'-k 

D*p 

Correlated background 

Uncorrelated background 



32844 ± 181 

21741 ±217 0.662 ±0.007 
2091 ± 45 0.064 ± 0.001 
2703 ± 42 0.082 ± 0.001 
6287 ±40 0.191 ±0.001 



D. At Fit Procedure 

In order to measure the CP violation parameters in 
the D*tt sample, we perform a simultaneous unbinned fit 
to the same-flavor (SF) events, in which the fast pion and 
the tagging lepton have the same charge, and opposite- 
flavor (OF) events, in which the fast pion and the tagging 
lepton have the opposite charge. We minimize the quan- 
tity — 21n£ = — 2J2t In A) where 

*~*i J D*tv-l D*7r T" JD*piD*p ± JuncoMinco "T Jcorr-fcorr- 

(11) 

The event-by-event signal and background fractions 
(the / terms) are taken from the results of the kinematic 
fits. Each P term contains an underlying physics PDF, 
with experimental effects taken into account. For D*ir 
and D*p, the PDF is given by Eq. 1, where for D*p the 
terms S are effective parameters averaged over the he- 
licity states [22] that are constrained to be zero. The 
PDF for correlated background contains a term for neu- 
tral B decays (given by Eq. 1 with S l± = 0), and, in the 
case of OF events, a term for charged B decays (for which 
the PDF is 75-^ — e -l A *l/ T B + where t b + is the lifetime of 

ZT B + 

the charged B meson). 

The PDF for uncorrelated background also contains 
neutral and charged B components, with the remainder 
from continuum e + e~ — > qq (q = u,d,s,c) processes. 
The continuum PDF is modelled with two components: 
one with negligible lifetime, and the other with a finite 
lifetime. The sideband parameters are determined from 
data sidebands, as described later. 

As mentioned above, experimental effects need to be 
taken into account to obtain the P terms of Eq. 11. 
Mistagging is taken into account using 

(1 - w-)P(B° -> D^w*) + w+P(B° -> D**^) 

(1 - w+)P(B° -> ZTTtt*) ± w_P(B° -> Zr^TT*), 

(12) 



where w+ and w_ are the wrong tag fractions, and are 
determined from the data as free parameters in the fit 
for 5 f± . It should be noted that the w+ and w- values 
used here are different from those used in the full recon- 
struction methods because the flavour-tagging methods 
differ in two cases. 

The time difference At is related to the measured quan- 
tity Az as described in Eq. 4, with an additional term 
due to possible offsets in the mean value of Az, 



At 



At + e At ~(Az + e Az )//3~fc. 



(13) 



It is essential to allow non-zero values of e since a small 
bias can mimic the effect of CP violation: 

cos(AmAt) -> cos(AmAt) - Ame A t sin(AmAt). (14) 

A bias as small as e Az ~ 1 /im can lead to sine-like terms 
as large as 0.01, comparable to the expected size of the 
CP violation effect. We allow separate offsets for each 
combination of fast pion and tagging lepton charges. We 
also apply a small correction to each measured vertex po- 
sition to correct for a known bias due to the relative mis- 
alignment of the SVD and CDC in SVD1 data. This cor- 
rection is dependent on the track charge, momentum and 
polar angle, measured in the laboratory frame. It is ob- 
tained by comparing the vertex positions calculated with 
the alignment constants used in the data, to those ob- 
tained with an improved set of alignment constants [24] . 
The alignment in SVD2 data was found to be compara- 
ble to that of the corrected SVD1 data. No additional 
correction was thus applied to SVD2 data. 



E. At Resolution 

Resolution effects are taken into account in a way sim- 
ilar to our full reconstruction analysis. The algorithm 
includes components related to detector resolution, kine- 
matic smearing and non-primary tracks. 

For correctly tagged signal events, both the fast pion 
and the tagging lepton originate directly from B me- 
son decays. Therefore we do not include any additional 
smearing due to non-primary tracks in these events. In- 
correctly tagged events, however, almost exclusively orig- 
inate from secondary leptons. Furthermore, due to the 
kinematic constraints on the momentum of the fast pion, 
secondary tracks consist almost exclusively of wrong-tag 
leptons. Only uncorrelated background events contain a 
small amount of secondary pions, which also give the 
wrong flavour information. In order to take this ef- 
fect into account, the PDFs of incorrectly tagged events 
are convolved with an additional resolution component 
whose parameters arc determined from MC simulations. 
Three different sets of parameters are used: one set for 
uncorrelated background in the uncorrelated background 
sideband, where the fast pion momentum cut is loos- 
ened; one set for uncorrelated background in the two 
other regions (signal region and correlated background 
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side-band) ; one set for all other categories of events (sig- 
nal, D* p and correlated background), which contain sim- 
ilar amounts of secondary leptons, and no secondary pi- 
ons. Because of the aforementioned correlation between 
mistagging and non-primary tracks, each of these sets are 
linked to different wrong-tag fractions w± . 

The effect of the approximation that the B mesons are 
at rest in the cms in Eq. 4 is taken into account [23] . We 
use a slightly modified algorithm to describe the detector 
resolution, in order to precisely describe the observed be- 
haviour for single track vertices. The resolution for each 
track is described by the sum of three Gaussian compo- 
nents, with a common mean of zero, and widths that are 
given by the measured vertex error for each track multi- 
plied by different scale factors. 

We measure the five parameters of the detector reso- 
lution function (three scale factors and two parameters 
giving the relative normalizations of the Gaussians) us- 
ing J/ip — > p + p~ candidates. These are selected using 
criteria similar to those for D*tt, except that both tracks 
are required to be identified as muons, and their invari- 
ant mass is required to be consistent with that of the 
J/tp- Vertex positions are obtained independently for 
each track, in the same way as described above. The 
quantity Az — z^+ — z^- then describes the detector 
resolution, which is the convolution of the two vertex 
resolutions, since for J/ip — > p + pT the underlying PDF 
is a delta function. 

We perform an unbinned maximum likelihood fit using 
events in the J/ip signal region in the di-muon invariant 
mass, and using sideband regions to determine the shape 
of the background under the peak. The underlying Az 
PDF of the background is parametrized in the same way 
as that used for continuum, described above. Two sets of 
parameters arc simultaneously determined for SVD1 and 
SVD2 data. We also take a possible offset into account, 
so that there are in total 16 free parameters in this fit 
(five describing the resolution function, two describing 
the background, and one Az offset, each for SVD1 and 
SVD2 data). The result is shown in Fig. 7. 



F. Background At Shape 

The free parameters in the background PDFs P U nco 
and P C orr are determined using data from sideband re- 
gions. To measure the uncorrelated background shape, 
we use events in a side-band region of —0.5 < cos 5 7Tf7rs < 
0.5, which by definition is only populated by uncorre- 
lated background, thus providing a very large sample to 
extract the various components of this background. We 
perform a simultaneous fit to OF and SF candidates, in 
both SVD1 and SVD2 data. 

To obtain the correlated background parameters, a si- 
multaneous fit is carried out to OF and SF events in a 
sideband region of 0.85 < cosS„,„ e < 0.90, 1.93 GeV/c < 
p Wf < 2.13 GeV/c and -1.7 < cos6» hc i < -0.3. This side- 
band region is dominated by correlated and uncorrelated 
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FIG. 7: Result of the resolution parameter extraction proce- 
dure for J/tp — > p + p~ candidates selected from (top) SVD1 
and (bottom) SVD2 data, shown with both (left) linear and 
(right) logarithmic coordinate scales. The data points show 
the Az distribution for candidates in the signal region; the 
full line shows the result of the fit. The dashed line indicates 
the background. 



backgrounds — the contributions from D*ir and D* p arc 
found to be small in MC. The uncorrelated background 
parameters are fixed to the values found in the previous 
fit, except the Az offset (which are the same as that of 
correlated background) and the wrong-tag fractions, as 
mentioned above. 



G. Fit Result 

In order to test our fit procedure, we first constrain 
S + and S~ to be zero and perform a fit in which tqo 
and Am (as well as two wrong tag fractions and eight 
offsets) are free parameters. We obtain tb« — 1.495 ± 
0.012 ps and Am = 0.506 ± 0.006 ps _1 , where the errors 
are statistical only. These values are compatible with the 
current world averages [19]. Reasonable agreement with 
the input values is also obtained in MC. Furthermore, fits 
to the MC with S floated give results consistent with 
zero, as expected. 

To extract the CP violation parameters we fix t b o 
and Am at their world average values, and fit with 
5 + , S~ , two wrong tag fractions and eight offsets as 
free parameters. We obtain S + = 0.048 ± 0.028 and 
S~ = 0.034 ± 0.027 where the errors are statistical only. 
The wrong tag fractions are w_ = (4.8 ± 1.6)% and 
w + = (3.4 ± 1.6)%. All floating offsets are consistent 
with zero except the one for the ir~£ + combinations in 
the SVD1 sample. The results are shown in Fig. 8. To 
further illustrate the CP violation effect, we define asym- 
metries in the same flavour events (»4 SF ) and in the op- 
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posite flavour events („4 OF ), as 



H. Systematic Error 



^ SF 
^ OF 



N v - l -(Az)-N„ +l+ (Az) 
N„- l -(Az)+N„ +l +(Az) 
N„ +l -(Az)-N„- l+ (Az) 

N v+l -(Az) + N^- l+ (Az) 



(15) 



where the N values indicate the number of events for each 
combination of fast pion and tag lepton charge. These 
are shown in Fig. 9. Note that due to the relative con- 
tributions of the sine terms in Eq. 1, vertex biases (i.e. 
non-zero offsets) can induce an opposite flavour asymme- 
try, whereas the same flavour asymmetry is more robust. 
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FIG. 8: Results of the fit to obtain 5+ and S~ . The fit result 
is superimposed on the data. The signal component is shown 
as the dashed line. The dotted line indicates the background 
contribution. 



Systematic errors due to the resolution function, the 
background fractions, the background parameters, and 
physics parameters are estimated by varying the values 
used in the fit by ±ler. Systematic effects due to differ- 
ences between data and MC in the distributions used in 
the kinematic fit are further investigated by repeating the 
fit using different binning. We repeat the entire fit pro- 
cedure using twice as many bins in each of the three dis- 
criminating variables. Since cos<5 T/ ^ s is used in the best 
candidate selection, we also repeat the algorithm with- 
out using this variable in the kinematic fit. The largest 
deviation (0.007) is assigned as a systematic error. 

The resolution function parameters are precisely de- 
termined from the fit to J/ip — > p + p~ candidates. We 
consider systematic effects due to our lack of knowledge 
of the exact functional form of the resolution function: 
using different parametrizations results in shifts of S 
as large as 0.008, which we assign as a systematic error. 
Allowing for effective 5 terms of ±0.05 in the D*p and 
correlated background PDFs leads to systematic errors 
of 0.003 and 0.002, respectively. 

Vertex biases may lead to a large systematic error on 
S ± , so we have introduced offsets to make our analysis 
relatively insensitive to such biases. These additional free 
parameters cause an increase in the statistical error of 
about 20%. In order to estimate the systematic error due 
to these offsets, we compare the results in MC simulations 
with and without the offsets. A difference of +0.011 is 
found for S + and is assigned as a systematic error for 
the offsets. We have further tested our fit routine for 
possible fit biases, such as could be caused by neglecting 
terms of 11%,^ in the fit, by generating a number of large 
samples of signal MC with different input values of S + 
and S~ . All results come out consistent with the input 
values, without evidence of any bias. 

The systematic errors are summarized in Table V. The 
total systematic error (0.017) is obtained by adding the 
above terms in quadrature. 
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FIG. 9: Results of the fit to obtain S and S~ , shown as 
asymmetries in the (left) same flavour events and (right) op- 
posite flavour events. The fit result is superimposed on the 
data. 



TABLE V: Summary of the systematic errors in the S mea- 
surements using the partial reconstruction method. 



Source 



Error 



Resolution fit 0.002 

Resolution models 0.008 

Kinematic smearing 0.002 

Non-primary tracks 0.004 
Background shapes negligible 

Kinematic fit 0.007 

T B o,Am 0.001 
CP violation in D* p and corr. bkgd. 0.004 

Vertexing 0.011 



Total 



0.017 
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I. Result 

The results using the partial reconstruction method 
are 

S + = 0.048 ±0.028 ±0.017, 

S~ = 0.034 ±0.027 ±0.017, (16) 

where the first error is statistical and the second error is 
systematic. 
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V. DISCUSSION 

A. Final results of S ± (D*tt) and S ± (Dtv) 

The D*ir signal candidates in the full reconstruction 
sample and the partial reconstruction sample are mostly 
independent. We find only 60 D*ir candidates which en- 
ter both samples. We repeat the analysis of the partial 
reconstruction sample after removing those overlapping 
candidates and obtain the same result. Therefore we 
combine the two results. Some part of the systematic 
errors in the two measurements may be correlated. We 
conservatively assume that contributions from physics 
parameters and fit biases (fit and MC bias in the full 
reconstruction sample and Az offsets in the partial re- 
construction sample) in the two measurements are 100% 
correlated, and combine them by taking a weighted aver- 
age using the inverse of the statistical errors as weights. 
The final results are 

S + (D*ir) = 0.049 ±0.020 ±0.011, 

S-(D*ir) = 0.031 ±0.019 ±0.011, 

S+(Dtt) = 0.031 ±0.030 ±0.012, 

S-(Dir) = 0.068 ±0.029 ±0.012 



(17) 



where the first errors are statistical and the second errors 
are systematic. These results are shown in Fig. 10 in 
terms oil a, 2 <r and 3 a allowed regions in the S~ versus 
S + space. Small possible correlations in the systematic 
errors of S + and S~ are neglected. 

The results for D*ir and Dir decay modes show that 
the (S + + S~) values are 2.5 a and 2.2 a away from 
zero while the (S + — S~) values are within 1 a of zero. 
Since (S + + S~) oc sin(20i + 3 )cosc5 and (5+ - S~) oc 
cos(20i + 03) sin 5 (Eq. 2), it can be seen that these re- 
sults are consistent with both Sd'tt and 5^^ being small, 
as predicted by some theoretical models [5]. The signif- 
icance of CP violation, seen as deviations of (5+ + S~ ) 
from zero, is 2.5 a for D*ir and 2.2 a for Dir decay modes. 



B. Constraints on (20i + 03) and R 

Since we have two measurements {S + and S~ ) which 
depend on three unknowns (i?, 20i + 03, 5), there is 



FIG. 10: Results of the S measurements expressed in terms 
of S~ vs S + for the D*tt (left) and Dir (right) modes. Shaded 
regions indicate allowed regions with la, 2a and 3 a uncer- 
tainties defined by y/— 21n£ = 1, 4, 9, respectively 



not sufficient information to solve for the weak phase 
(20i +03). Instead we obtain exclusion regions in two- 
dimensional space for any value of the third variable. The 
regions of (20i + 03, R) space that are excluded at the 
1, 2, 3cr levels are shown in Fig. 11. An alternative 
representation, shown in Fig. 12, gives the lower bound 
on J sin(20i + 03 )| for any values of R and 6. 
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FIG. 11: Excluded regions of(20i +03 ) vs R space at 1, 2 , 3 a 
level for the D*n (left) and Dir (right) decays. The range 0°- 
180° for 20i + 03 is shown; there are additional solutions at 

201+03 — > 201 +03 + 180°. 

Further conclusions cannot be drawn without some 
theoretical estimate of the values of either R D ^)^ or 
8 D {,)^. One interesting possibility is to estimate the 



of R 



DM-, 



using decays such as i? H 



B {) 



D s n T , which are related to B' 



D^^tt and 
* ,D(*)±7rF by 
isospin and SU(3), respectively [2]. 

The method using SU(3) symmetry has some experi- 
mental advantages, since the rates are enhanced by the 
square of the tangent of the Cabibbo angle 9 C . The rele- 
vant expression is 



R 



D(*)tt 



= tan 



/. 



D" 



B(B° ^ Dl* ,+ ir-) 



f 



-D 



B(B° -> L>(* 



")' 



(18) 



where Jm is the decay constant for the M -meson. The 
equality is valid up to SU(3) breaking effects. Both 
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FIG. 12: Lower limit on | sin(2c6i + <fe)| as a lunction of R at 
lcr (68% CL) (thick curves) and 2a (95% CL) (thin curves) 
from (left) D*n and (right) Dn. The estimated values for 
R D (t) n are indicated as solid lines (central values). Coarse 
dotted lines and fine dotted lines indicate their lcr and 2a 



Belle [26] and BaBar [27] 
the decay B° -> D+tt". 



have reported evidence for 
BaBar included a limit for 



the B° — ► D* + tt~ mode in the same paper. We es- 
timate Rd*tt and Rd-k using the measured branching 
fractions and the form factors from a lattice QCD cal- 
culation by the UKQCD collaboration [28]. We use 



B(B° -> D* s + n-)/B(B° 
which is determined from 



D* 



0.0068] 



-0.0051' 

the BaBar result B(B° -> 
, -TT-) ---- (1.9±i| ± 0.5) x 10~ 5 and the PDG 2004 
value B(B° -> D*~-k + ) = (2.76 ± 0.21) x 10~ 3 , and 
B(B° -> Dfn-)/B(B -> D-JT+) = 0.0098 ± 0.0040, 
which is determined from the PDG 2004 values of 
B(B° -> D+n-) = (2.7 ± 1.0) x 10~ 5 and B(B° -> 
D--K+) = (2.76 ± 0.25) x 10~ 3 . The decay constant ra- 
tios are JW/d* = 1-04 ± 0.01 ± 0.02 and f D /f D . = 
0.90 ±0.01 ±0.01. Weusetan6» c = 0.23 [19, 29]. 

Adding ±30% theory uncertainty for SU(3) breaking 
effects, such as contributions from annihilation diagrams, 
we obtain 

R D . T = 0.020 ± 0.007 ±0.006(theory), 
R D „ = 0.021 ± 0.004 ±0.006(theory) (19) 

as indicated in Fig. 12. If we use these values, we can 
restrict the allowed region of | sin(20i + 03) | to be 

| sin(20i + 03)| > 0.44 (0.13) at 68%(95%) CL (20) 
from the results for the D*tt mode and 

| sin(20i + 3 )| > 0.52 (0.07) at 68%(95%) CL (21) 
from the results for the Dtt mode, respectively. 

VI. SUMMARY 

We have measured CP violation parameters that de- 
pend on 03 using the time-dependent decay rates of the 



B° -> flW^. A total of 386 million BB events were 
used in the analysis. While the Dtt sample was collected 
by a standard full reconstruction method, the D*ir sam- 
ple was enlarged by using, in addition, a partial recon- 
struction technique, essentially doubling the statistical 
power of this mode. 

The final results expressed in terms of S + and S~, 
which are related to the CKM angles 0i and 03, 
the ratio of suppressed to favoured amplitudes, and 
the strong phase difference between them, as 5 ± = 
-R D '« sin(20i + 03 ± <W)/ (1 + Rd**) for D *^ and 
S ± = +Rd-k sin(20i + 03 ± 6 D „)/ (l + R%„) for Dtt, are 



S + (D*n) = 0.049 ±0.020 ±0.011, 

S-(D*tt) = 0.031 ±0.019 ±0.011, 

S + (Dir) = 0.031 ±0.030 ±0.012, 

S-(Dtt) = 0.068 ±0.029 ±0.012, 



(22) 



where the first errors are statistical and the second errors 
are systematic. These results are an indication of CP vi- 
olation in B° — > D*~tt + and B° — > D~ir + decays at the 
2.5 a and 2.2 c levels, respectively. If we use the values of 
Rd'-k and Rd-k that are determined using a combination 
of factorization and SU(3) symmetry assumptions, the 

branching fraction measurements for the D s it modes, 
and lattice QCD calculations, we obtain 68% (95%) con- 
fidence level lower limits on | sin(20i + 03)| of 0.44 (0.13) 
and 0.52 (0.07) from the D*ir and Dtt modes, respec- 
tively. 
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